Pattern formation of quantum jumps with Rydberg atoms 
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We study the nonequilibrium dynamics of quantum jumps in a one-dimensional diain of atoms. 
Each atom is driven on a strong transition to a short-hved state and on a weak transition to a 
metastable state. We choose the metastable state to be a Rydberg state so that when an atom 
jumps to the Rydberg state, it inhibits or enhances jumps in the neighboring atoms. This leads to 
rich spatiotemporal dynamics that are visible in the fluorescence of the strong transition. It also 
allows one to dissipatively prepare Rydberg crystals. 
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Introduction. Pattern formation is a universal phe- 
nomenon in nonequilibrium systems [ij. In the presence 
of driving and dissipation, a system may spontaneously 
form sustained spatial structures. Pattern formation is 
usually studied in classical systems, such as fluids, chem- 
ical reactions, and granular media. But can a many-body 
quantum system exhibit pattern formation? This is not 
obvious, because adding dissipation to a quantum system 
usually leads to decoherence, so one might not expect any 
interesting behavior. Also, quantum mechanics is linear, 
while pattern formation in classical systems requires non- 
linearity. 

In this paper, we show that a quantum system of Ryd- 
berg atoms can form spatial patterns in its fluorescence. 
A Rydberg atom is an atom excited to a high princi- 
pal quantum number n, and there is a strong dipole- 
dipole interaction between nearby Rydberg atoms. This 
interaction leads to a wide range ofphysics, including 
quantum information processing [2|-l5|, quantum phase 
transitions mp, thermalization of closed quantum sys- 
tems nonlinear optics (l^ - [T5l |. and dissipative 
physics pii, ,17]. Here, we show that Rydberg atoms also 
support pattern formation. 

Our idea is based on quantum jumps of a three-level 
atom p^| - |20j . It is well known that an atom driven 
strongly to a short-lived state and weakly to a metastable 
state occasionally jumps to and from the metastable 
state. The jumps are visible in the fluorescence of the 
strong transition, which exhibits distinct bright and dark 
periods [2T1 - [23| . Quantum jumps are a good example of 
how a quantum system far from equilibrium (due to laser 
driving and spontaneous emission) can have nontrivial 
dynamics. 

Here, we consider a lattice of many three-level atoms, 
and we let the metastable state be a Rydberg state, so 
that a jump of one atom affects its neighbors' jumps 
via the dipole-dipole interaction. This leads to rich 
spatiotemporal dynamics in a one-dimensional chain of 
atoms. Depending on the experimental parameters, there 
is a variety of behaviors. The patterns are observable by 
imaging the fluorescence of the strong transition. Our 
scheme can also be used to make Rydberg crystals [B-Q ■ 

Previous works studied correlated quanturn jumps of 
atoms in the context of the Dicke model [13, ■ They 
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FIG. 1: (a) An atom has a ground state \g), short-lived ex- 
cited state |e), and metastable state |r), which is chosen to 
be a Rydberg state. One observes the spontaneous emission 
from |e). (b) The \g) o \r) transition is originally on reso- 
nance (Ar = 0), but when one atom is in \r), the other atom 
is off resonance, (c) The \g) -f^ \r) transition is originally off 
resonance (Ar = V), but when one atom is in |r), the other 
atom is on resonance, (d) When Ar = 0, \rr) is weakly cou- 
pled to the other states. Note that (b) and (d) are equivalent. 



concluded that cooperative effects are very difficult to see 
experimentally, because the interatomic distance must be 
much smaller than a wavelength. In contrast, the strong 
Rydberg interaction here allows the interatomic distance 
to be much longer than a wavelength. Thus, the atoms 
develop strong correlations while being individually re- 
solvable. 

Results for a single atom. We first review quantum 
jumps in a single atom [18. .20,] . Consider an atom with 
three levels: ground state \g), short-lived excited state 

e), and metastable state \r) [Fig. dja)]. We choose the 
metastable state to be a Rydberg state since Rydberg 
states have long lifetimes [2g[. A laser drives the strong 
transition \g) ^ |e), while another drives the weak tran- 
sition \g) o \r). The strong transition acts as a mea- 
surement of whether or not the atom is in |r). When 
the atom is not in |r), the atom is repeatedly excited to 

e) and spontaneously emits photons. Occasionally the 
atom is excited to |r) and stays there, and the fluores- 
cence from the strong transition turns off. Eventually, 
the atom returns to and the fluoresence turns back 
on. Thus, the fluorescence signal of the strong transition 
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exhibits bright and dark periods, and the occurrence of 
a dark period implies that the atom is in \r). 
The Hamiltonian for a single atom is {fi= 1) 



H = 



-Ae|e)(e| 



\e){9\) 
Ar\r){r\ 



{\9){r 



0(.9l) 



(1) 



where Ag and fie are the laser detuning and Rabi fre- 
quency of the strong transition, while and fir are 
the corresponding quantities for the weak transition. In 
the absence of spontaneous emission, Eq. ([!]) would com- 
pletely describe the system. However, the excited states 
have lifetimes given by their linewidths, 7e and jr- 

In the rest of paper, we make the following assumptions 
on the parameters. To avoid power-broadening on the 
strong transition, we choose to work in the low-intensity 
limit, He ^ 7e; this choice is clarified later. For conve- 
nience, we set Ae = 0, although it may be experimentally 
useful to set Ae < for continuous laser cooling [23| . We 
also set 7r. = 0, since the lifetime of the Rydberg state 
scales as and hence can be chosen to be arbitrarily 
long [2^. It is straightforward to extend the analysis to 
nonzero Ae and 7^. 

Well-defined jumps appear in the fluorescence signal 
when a bright period consists of many photons while a 
dark period consists of the absence of many photons. For 
a single atom, this happens when ftr ^ ^llle in the case 
of Ar = [19[. The transition rate from a dark period 
to a bright period is (20j 



(A.) = 



16A4-|-4A2(72_2f72) + f7;! 



(2) 



and the rate from a bright period to a dark period is 
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where B and D denote bright and dark periods. Both 
rates are maximum when A^ — since the strength of 
the weak transition is maximum there. When A^ = 0, 
both rates are approximately 7eO^/rig. This depends 
inversely on fig, because increasing fig is equivalent to 
measuring the atomic state more frequently; this inhibits 
transitions to and from |r), similar to the quantum Zeno 
effect m. 

Many-body model. Now we consider a one-dimensional 
chain of N three-level atoms. The atoms are coupled 
via the dipole-dipole interaction between their Rydberg 
states. In the absence of a static electric field, the inter- 
action decays with the sixth power of distance for long 
distances [5| ; hence, we assume nearest-neighbor interac- 
tion. The Hamiltonian is 
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FIG. 2: Quantum trajectory simulation of chain of A'' = 8 
atoms with periodic boundary conditions. The Rydberg pop- 
ulation of each atom is plotted vs. time, using color scheme 
on the right. White color means that the atom is bright and 
not in the Rydberg state. Black color means that the atom is 
dark and in the Rydberg state, (a) ~ 0.2^e, Q,- = 0.0057e, 

A,. = 0, V ^ 0.l7e. (b) fie = 0.l7e, = 0.0057e, 

Ar = V = 0.l7e. (c) = = 0.l7e, A,. = 0, = 0.37e. 



where V is the nearest-neighbor interaction and we as- 
sume periodic boundary conditions. The interatomic dis- 
tance is assumed to be large enough so that the fluores- 
cence from each atom is resolvable in situ on a camera 

To demonstrate the rich spatiotemporal dynamics of 
the many-body system. Fig. [5] shows simulations of a 
chain of = 8 atoms, generated using the method of 



quantum trajectories ^ 3lj- Each trajectory simulates 



a single experimental run. We have plotted the Rydberg 
population of each atom, i.e., the expectation value of 
Ri = \r){r\i. The atoms undergo quantum jumps, and 
the Rydberg interaction clearly leads to spatial correla- 
tions in the fluorescence. 

There are different types of collective dynamics 
depending on the parameters. In Fig. [2ja)-(b), 
fir <C rig/7e, so each atom by itself would exhibit quan- 
tum jumps. In Fig. [IJa) (A^ = 0), a dark period usually 
does not spread to the neighboring atoms. But once in 



3 



0.25 




0.75 



FIG. 3: Quantum trajectory simulation of A*' = 2 atoms. The 
Rydberg population of each atom is plotted vs. time, using 
color scheme on the right. White color means that the atom is 
bright and not in the Rydberg state. Black color means that 
the atom is dark and in the Rydberg state. Parameters are 
the same as in Fig. (2] (a) = 0.27^, Qr = 0.0057e, Ar = 0, 
V = 0.l7e. (b) f7e = 0.27e, Q.r = 0.0057e, = y = 0.l7e. 

(c) ile^^r= 0.l7e, A^ = 0, F = 0.37e. 



a while, there are two dark atoms in a row (BDDB) 
and they both stay dark for a long time. In Figs. [UJb) 
(Ar = V), there are two types of behavior. Sometimes, 
dark periods spread across the chain so that many con- 
secutive atoms are dark at the same time. Other times, 
an alternating DBDB pattern appears. Over time, the 
system switches between these two types of behaviors. 
In Fig. ^c) {fir = r^e, A^ = 0), the atoms tend to turn 
dark or bright in groups of two or three, and sometimes 
all the atoms are dark. The existence of jumps here is 
surprising because a single atom would not exhibit jumps 
for these parameters. 

To understand the results for = 8, it is instructive 
to consider the simpler case of = 2 atoms. Figure [3] 
shows quantum trajectory simulations for TV = 2; note 
the similarity with Fig. [5] We have analytically solved 
the N = 2 case. The details are in the supplemental ma- 
terial [111 , and here we summarize the results and provide 
intuition. There are two general cases: (i) ftr ^ ^l/le 
and (ii) 51^ = J7e, A^ = 0, distinguished by whether a 
single atom would exhibit jumps. 

Case of fir ^ r2g/7e. For these parameters, an atom 
by itself would exhibit jumps. Let the two atoms be 
labelled 1 and 2. If atom 1 is in |r), then according to 
Eq. dH, atom 2 effectively sees a laser detuning of Ar — V. 
But if atom 1 is not in |r), then atom 2 sees the origi- 
nal detuning A^. Whether atom 1 is in |r) depends on 
whether it is in a dark period. This suggests that the 
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FIG. 4: Statistics for TV = 2 atoms, (a) Ratio of r""^"" to 
yBD^BB j-j^ ^ Q 27^, = 0.0057e, V = 0.l7e. (b) Jump 
rates for ilr = Qe = 0.l7e, A^ = 0. T^o^bb. analytical 
result (black, solid line) and numerical data (black circles). 
Y^BB^DD, analytical upper bound (blue, dashed line) and 
numerical data (blue triangles). 



jump rates for atom 2 are the same as for a single atom 
[Eqs. ©-([S])], except with an effective detuning that de- 
pends on whether atom 1 is in a bright or dark period at 
the moment. In the supplemental material [12], we use a 
more careful analysis to show that this is indeed correct 
in the limit of small fir. Thus, the transition rates for 
two atoms are 
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An insightful quantity is the ratio r^^^-^^/r^^^"^^, 
which indicates how often the atoms are in DD relative 
to BB. As shown in Fig. |4Ka), the ratio is minimum at 
Ar = and maximum at A^ = V. 

The minimum at A^ = is due to the blockade effect: 
when atom 1 is in |r), it shifts the Rydberg level of atom 2 
off resonance so that atom 2 is prevented from jumping to 
k) [Fig-IHb)]. Thus, the atoms switch between SB, 
and DB; they are almost never in DD. In other words, 
there is at most one dark atom at a time [Fig. [2^a)]. 

The maximum at A^ = 1^ is due to the opposite ef- 
fect: the laser is originally off-resonant, but when atom 
1 happens to jump to |r), it brings the Rydberg level of 
atom 2 on resonance, encouraging atom 2 to jump to |r) 
[Fig. [TJc)]. Thus, the atoms switch between DD, BD, 
and DB; they are almost never in BB, except for the 
initial transient. When the atoms are in DD, there is an 
equal chance to go to BD or DB. Thus, the dark spot 
appears to do a random walk between the two atoms 
[Fig.Hb)]. 

The above considerations also apply to larger TV. In 
general, the transition rates for atom i are given by 
Eqs. (H])-© but with an effective detuning that depends 
on its neighbors' current state: A^ — V{{Ri-i) + {Ri+i)). 
When Ar — 0, the blockade effect prevents dark peri- 
ods from spreading [Fig. [UJa)]. However, if there hap- 
pens to be two dark atoms in a row (BDDB), the dark 
atoms are effectively off resonance and are unlikely to 
turn bright; thus, they stay dark for a long time. On 
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the other hand, when = V, the anti-blockade ef- 
fect encourages dark periods to spread across the chain 
[Fig. [2ljb)] . However, sometimes an alternating pattern 
develops (DBDB); this is because when an atom is 
bright while both neighbors are dark, it is effectively off 
resonance and is unlikely to turn dark. 

Case of fir = ^e, — 0. For these parameters, an 
atom by itself would not exhibit jumps because of the 
absence of a weak transition. The existence of jumps for 
two atoms is soley due to the dipole-dipole interaction, 
which causes \gr) o \rr) and \rg) 4t- \rr) to become 
off-resonant and thus weak transitions [Fig. [TJd)]. Since 
\rr) is metastable, occasionally the system jumps to and 
from \rr). When the system is in |rr), the atoms do not 
fluoresce. When the system is not in |rr), the wavefunc- 
tion rapidly oscillates among the other eigenstates so that 
both atoms fluoresce from |e). Thus, the system switches 
between BB and DD [Fig. [Sljc)]. The rates are ^ 
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YBB^DD ^ 
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where = fi^ = 17e. The inequality for r^^^-°^ is due 
to incomplete knowledge of the wave function after a pho- 
ton emission. Equations (|9|- (fT0| agree well with quan- 
tum trajectory simulations [Fig. IHb)]. Both rates are 
inversely related to V , since the weak transitions become 
weaker as V increases. The condition for well-defined 
jumps is roughly ^ <^2V. 

A larger chain seems to have similar behavior 
[Fig. IHc)]. The atoms tend to turn dark or bright simul- 
taneously with their neighbors. However, the dynamics 
are more complex due to the presence of two neighbors. 

Experimental considerations. Our scheme can be im- 
plemented using atoms trapped in an optical lattice. For 
example, one can use ®''Rb, which has a strong hS — 5P 
transition with linewidth 7e/27r = 6 MHz [27i] . Suppose 
one chooses the 605 Rydberg state, which can be reached 
via a two-photon transition. For a lattice spacing of 



7 /xm, the nearest- neighbor interaction is 1/ = 0.27e 
The lifetime of that Rydberg state is 250 /j,s at K 
in other words, 7^ « 7e/10'*. Transitions due to black- 
body radiation can be minimized by working at cryogenic 
temperatures. Also, the nS states have neg ligible losses 
from trap-induced photoionization [ssl l36j | The trap- 
ping of Rydberg atoms in optical lattices was recently 
demonstrated in Refs. [36, 37]. 

There is an important constraint on the experimen- 
tal parameters: the interaction V should be much less 
than the trap depth, or else the repulsive interaction will 
push the atoms out of the lattice. Since a trap depth 
of 10 MHz is possible [slj, we require F <C 7e- Then to 
avoid broadening the strong transition [l9j and smearing 
out the effect of V, we choose fie ^ 7e, as stated above. 

Quantum jumps can also be used to prepare a Ryd- 
berg crystal, which is a regular arrangement of Rydberg 
atoms among ground-state atoms (6|-|9|. Rydberg crys- 
tals are important for studying many-body physics and 
generating non-classical light, but are nontrivial to make. 
Here, one would let the atoms jump to and from the Ry- 
dberg state while monitoring the fluorescence until the 
desired crystal is obtained. The parameters can be opti- 
mized for the desired crystal type. This is in the spirit 
of recent works that use dissipation to prepare nontrivial 
quantum states [H, [s^ 

Conclusion. Thus, quantum jumps of Rydberg atoms 
lead to interesting spatioteniporal dynamics. The next 
step is to see what happens in larger systems, especially 
in higher dimensions: does the Rydberg interaction lead 
to long-range patterns in a large system? It would also 
be interesting to see what happens with other interaction 
types (long range or anisotropic) or when the atoms are 
free to move instead of being fixed on a lattice. Finally, 
one should look for other quantum many-body systems 
that exhibit pattern formation. Perhaps nonequilibrium 
quantum systems share some universal features as is the 
case with nonequilibrium classical systems 

We thank H. Haffner and H. Weimer for useful dis- 
cussions. This work was supported by NSF Grant No. 
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photon by time given that it emitted at time 0. Po{t) 
decreases monotonically as t increases. When the param- 
eters are such that there are well-defined quantum jumps, 
Po(t) decreases rapidly to a small value for small i, but 
has a long tail for large t. This reflects the fact that the 
time between emissions is usually short (bright period), 
but once in a while it is very long (dark period). Note 
that each emission is an independent event, due to the 
fact that the wave function always returns to \g) after an 
emission. 

We write Po{t) = -Pshoit(i) + P\ong{t) to separate the 
short and long time-scale parts. The long tail is given 
by P\ong,{t) = pexp(— r^~*'^t), where p is the probability 
that a given interval is long enough to be a dark period, 
and r^^^ is the transition rate from a dark period to 
a bright period. In other words, is the average 

duration of a dark period. 

To calculate Pq {t) , we follow the evolution of the wave 
function |^/'(i)), given that the atom has not emitted a 
photon yet. This is found by evolving \ip{t)) with a non- 
Hermitian Hamiltonian Hog = H — i^\e){e\. The non- 
Hermitian term accounts for the population thats emits 
a photon, hence dropping out of consideration 1]. Thus, 

Poit) = m)\m)- 

In the basis {\g), |e), jr)}, the matrix form of Hcs is 



SUPPLEMENTAL MATERIAL 

This supplement provides details of the jump-rate 
calculations. Appendix A reviews the derivation for 
one atom. We essentially reproduce the derivation in 
Refs. [IHll, because we need to refer back to it later, and 
it is convenient to see it in our notation. Appendix B 
and C derive the results for two atoms: Appendix B cov- 
ers the case of ^Ir ^ ^e/7e; while Appendix C covers 
the case of fi^ = ^e, — 0. In general, we use the 
"quantum-trajectory" approach, which is based on the 
wave function, to account for spontaneous emission. An 
alternative approach is to use the density matrix as in 
Refs. I,!]. 



Appendix A: Review of one-atom case 

When an atom exhibits quantum jumps, the fluores- 
cence signal has bright periods, in which the photons are 
closely spaced in time, and dark periods, in which no pho- 
tons are emitted for a while. The goal is to calculate the 
transition rate from a bright period to a dark period and 
vise versa. The important quantity is the time interval 
between successive emissions [l[. During a bright period, 
the intervals are short, but a dark period is an exception- 
ally long interval. Suppose one has the function Po(^)j 
which is the probability that the atom has not emitted a 
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We want to solve the differential equation i-^\^p{t)) — 
Hcff\ip{t)) given the initial condition |V'(0)) = \g)- The 
general solution is \ip{t)) — X^n c„e~*'*'"*|u„), where 
and \un) are the eigenvalues and eigenvectors of -ffeff, 
and c„ are determined from the initial condition \g) = 

E„ C„\Un). 

We calculate the eigenvalues and eigenvectors pertu- 
batively in fir, which is assumed to be small. (Note that 
since Hcs is non-Hermitian, perturbation theory is differ- 
ent from the usual Hermitian case _6].) All three eigen- 
values have negative imaginary parts, which leads to the 
nonunitary decay. It turns out that the imaginary part 
of one of the eigenvalues, which we call A3, is much less 
negative than the other two. This means that the \ui) 
and \u2) components in \tpit)) decay much faster than 
the jus) component. After a long time without a photon 
emission, li^it)) contains only jua). Thus, A3 corresponds 
to the long tail of Po{t). 

To second order in il^ [E 13 1 



n^(-2A^-fz7e) 
8A2 - 2ni - 4i-feAr ■ 
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To first order in fir 
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Since [ua) consists mainly of |r}, tlie occurrence of a 
dark period implies, as expected, that the atom is in \r). 
(However, note that the atom is not completely in \r). In 
fact, the dark period ends when the small |e) component 
in |m3) decays and emits a photon Q.) 
We can now construct Piong{t): 
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Then instead of finding Pshort (t) explicity, we use a short 
cut 0. During a bright period, there is negligible pop- 
ulation in |r), so the atom is basically a two- level atom 
driven by a laser with Rabi frequency VLe- Thus, to low- 
est order in 12^, the emission rate Fghort during a bright 
period is the same as a two-level atom 



- short 



2m' 
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However, each emission in a bright period has a small 
probability p of taking a long time, in which case the 
bright period ends. Thus, the transition rate from a 
bright period to a dark period is 



P Tshort 
7e^+4A2 
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(AlO) 
(All) 



The jumps are well-defined when a bright or dark pe- 
riod is much longer than the typical emission time during 
a bright period: y^^^,T°^^ < Tshort- When A^ = 
and ^ 7e, this condition becomes 12^ ^ ^l/le- 



Appendix B: Two atoms, Q,r <^ ^i/le 

When Vlr <C $72/75, a single atom would exhibit quan- 
tum jumps. In the case of two interacting atoms, each 
one still undergoes quantum jumps, but the jump rates 
of each depend on the current state of the other atom. 
The goal is to calculate, to lowest order in 17^, the transi- 
tion rates among the possible states: BB^ BD, DB, and 
DD. 



Suppose for a moment that the interaction strength 
= 0. Then each atom jumps independently, and 
the jump rates are the same as the single-atom case 
[Eqs. ^Mi and (|IlT|) ]. 

Then let 7^ 0. Due to its form, the Rydberg inter- 
action only affects the state \rr). When the atoms are 
in BB, BD, and DB, there is negligible population in 
rr), so the interaction has negligible effect on the tran- 
sitions among BB, BD, and DB. So to lowest order in 
f7r, those transition rates are the same as when V ~ 0. 
Thus, we can immediately write down: 

YBB^bd _ Y^B^DB _ yB^d (J^^) 
YBD^bb _ Y^B^BB _ yD^b (Ji2) 

The remaining task is to calculate the transition rates 
that involve DD: T^d^^d _ yDB^dd _ y^d^bd ^ 
YDD^db 

To calculate these rates, we use an approach similar 
to Appendix A. Suppose the atoms are initially in BD, 
i.e., atom 1 is fluorescing while atom 2 is not. We are 
interested in the time interval between an emission by 
atom 1 and a subsequent emission by either atom 1 or 2. 
Usually the intervals are short since atom 1 is in a bright 
period. But once in a while, there is a very long interval, 
which means that atom 1 has become dark and the atoms 
are in DD. If the long interval ends due to an emission 
by atom 1, the atoms end up in BD; if it is due to an 
emission by atom 2, the atoms end up in DB. We want 
to calculate Poit), which is the probability that neither 
atom has emitted a photon by time t, given that atom 
1 emitted at time and also given that atom 2 started 
dark. Po(i) has a long tail corresponding to time spent 
in DD. 

We write Po{t) — -Pshoit(0 + Piong(i) to separate the 
short and long time-scale parts. The long tail is given 
by Piong(t) = pexp(-2r^^^-^^t), where p is the prob- 
ability that a given interval is long enough to be a DD 
period. 2T^^^^^ is the total transition rate out of DD 
since r^D^BD ^Y^D^DB. 

To evolve the wave function in the absence of an 
emission, we use the non-Hermitian Hamiltonian Hcs — 
H — i^(|e)(e|i + |e)(e|2), where H is the two-atom 
Hamiltonian. We want to solve the differential equa- 
tion i^lV'(O) = i?cff|V'(0) in order to find Po{t) = 

imim)- 

The question now is what initial condition to use. Since 
atom 1 is assumed to emit at time 0, it is in \g). Also, 
as discussed above, during a BD period, there is very 
little population in |rr), so the interaction has negligible 
effect on the dynamics. To first order in 17^, atom 2's 
wave function is the same as that of a single atom in a 
dark period [Eq. (IA3[) ]. So the initial condition of the 
two-atom system is: 



f7r(-2Ar -I- «7e) 



4A2 



172 _ 2i^^A, 



■\99) 



4A2 - 172 



2i7eA, 



-\9e) + \gr). (B3) 
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The general solution to the differential equation is 
[■^(O) = Z]nC„e"'^"*|M„), where A„ and |m„) are the 
eigenvalues and eigenvectors of Hcs, which is a 9 x 9 
matrix. c„ are determined from the initial condition 

We calculate the eigenvalues and eigenvectors pertur- 
batively in fir- All nine eigenvalues have negative imagi- 
nary parts, which leads to the nonunitary decay. It turns 
out that the imaginary part of one of the eigenvalues, 
which we call Ag, is much less negative than the other 
eight. This means that the other eight components of 
\ip{t)) decay much faster than the \ug) component. After 
a long time without a photon emission, lipit)) contains 
only jug). Thus, Ag corresponds to the long tail of Po{t). 

To second order in Qr, 



Ag 



-2Ar + V + 



172(_2A; + i7e) 



4a: 



J. (B4) 



where AJ, = A,. — V . To first order in il^, 
f7^(-2A; + z7e) 



|ug) 



4A;2 -ni- 2i7eA: 



4A;2 - 








-2A;h 




4A;2 - 













-|er) 



y\rg) 



cg 



4A;2 -ni- 2i7eA; 

»^(-2A; + ^7e) 

4A;2 -ni- 2z7eA; ' 



(B6) 



Note that |iig) consists mainly of |rr), since it corresponds 
to a DD period. 

We can now construct Piong(i): 



P 



|C9| 



16A;4 + 4A;2(72_2r!2)+174 

j^DD^BD ^ ^DD^DB = _ Ag 



16A;* + 4A;2(72_2172) + f74 



(B7) 
(B8) 
(B9) 
.(BIO) 



To calculate p-BC^UD^ ^-^^ short cut from Ap- 

pendix A. Since atom 1 is bright, it has negligible popu- 
lation in |r), so its emission rate Fghort is the same as a 
two-level atom [Eq. (IA9|) ]. Each emission has probability 
p of being long enough to be a dark period. 

j^BD^DD ^ r^S-^^=pr,hort, (Bll) 



7g + 4A;^ 



^DD^BD 



(B12) 



Note the similarity between Eqs. (IB10[) and (IA8|) and 
between Eqs. (IBT2)) and (UTTj) 



Appendix C: Two atoms, fir — Qe, A,. — 

For these parameters, a single atom would not exhibit 
quantum jumps. The existence of jumps for two atoms 
is solely due to the interaction. To calculate the jump 
rates, we use an approach similar to Appendices A and 
B, but there are some important differences. 

We are interested in the time intervals between photon 
emissions of either atom. We want to calculate Poit), 
which is the probability that neither atom has emitted a 
photon by time t, given that atom 1 emitted at time 0. 
(Alternatively, one could let atom 2 emit at time 0.) We 
write Po{t) = Pshort{t) + Piong{t) to separate the short 
and long time-scale parts. As in Appendix B, we want 
to solve the differential equation i-^\%p{t)) — i?off|V'(0) 
order to find Po{t) = {ip{t)\il){t)) . 

Before discussing what initial condition to use, wc first 
calculate the eigenvalues A„ and eigenvectors of iJcff- 
We define 17 = fi^ = and do perturbation theory in il, 
which is assumed to be small. As in Appendix B, all nine 
eigenvalues have negative imaginary parts, which leads to 
the nonunitary decay. The imaginary part of one of the 
eigenvalues, which we call Ag, is much less negative than 
the other eight. This means that the other eight compo- 
nents of I?/;) decay much faster than the |ug) component. 
Thus, Ag corresponds to the long tail of Po(i)- To fourth 
order in fi. 



Ag = y 
To first order in il, 

lug) 



2V 



■W) + TTTrlrg) + \rr), 



(CI) 



(C2) 



2V''' ' 2V' 

which consists mainly of |rr), refiecting the fact that if 
both atoms have not emitted for a while, they are in a 
DD period. 

Now it turns out that the real parts of the other eight 
eigenvalues have very different values, which causes the 
wave function to oscillate rapidly among the eight eigen- 
vectors. Thus, after atom 1 emits a photon, the short 
time-scale behavior consists of rapid oscillation among 
the eight eigenvectors, and each atom's |e) population 
fluctuates a lot. The time scale of the oscillation is faster 
than the typical photon emission rate, so both atoms are 
equally likely to emit next. Thus, the atoms can either be 
in BB or DD. When in BB, both atoms emit, and the 
time interval between emissions is relatively short. But 
once in a while, it takes a very long time for the next 
photon to be emitted, which means that the atoms are 
in DD. Once the long interval ends, the atoms go back 
to BB. 

The rapid oscillation during BB makes it impossible 
to choose a unique initial condition |?/'(0)), because each 
time atom 1 emits, atom 2's wave function is different. To 
account for this ignorance, we let atom 2's wave function 
be completely arbitrary: 



IV'(O)) = ai\gg) + a2\ge) + a3.\gr) . 



(C3) 
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Normalization requires |aip + |a2p + jasp = 1, but 
flij 02703 are otherwise unknown. Despite the incom- 
plete knowledge, we can still obtain a useful bound on 

YBB^dd 

The general solution to the differential equation 
i^|^/'(t)} = ifefflV-W) is IV^i)) = E„c„e-*^">„), where 
Cn are determined from the initial condition \ip{0)) = 
J2n Cn\un)- To first Order in il, 



as 



2V' 



(C4) 



Given the above results, we can now construct 
-Piong(^) = pexTp^—T^^^^^t), where p is the probability 
that a given interval is long enough to be a DD period, 
and r^^^-^^ is the transition rate from DD to BB: 



P 

j^DD-tBB 



|C9P < 



41/2 
-2 Im Ag 

21/2(^2+4^2)- 



(C5) 
(C6) 
(C7) 



The inequality for p reflects the incomplete knowledge of 



the initial wave function. 

To calculate p^s^dd^ have to first calculate Fshort, 
which is the total emission rate of both atoms during a 
BB period. We approximate Fshort using the emission 
rate in the absence of the \g) O \r) transition, like in 
Eq. jMI: 



short 



27e»' 

72 + 21^2 



(C8) 



However, since the \g) -H- |r) transition is not weak, the 
above approximation to Fshort is usually an upper bound. 
Now we can calculate: 



t — P t short 



< 



fi4 



(C9) 



The jumps are well-defined when a BB period consists 
of many emissions while a DD period consists of the ab- 
sence of many emissions: r^^^^^, F^^^^^ < Fshort- 
Roughly speaking, this happens when 



n <c 2V. 



(CIO) 
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